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ABSTRACT 

The  aluminum  oxide  crystal,  Al90~,  doped  with  a  small  amount  of 
chromium,  which  is  called  ruby,  is  widely  used  in  solid  state  physics 
research  and  applications.   The  experimental  work  described  in  this 
thesis  is  concerned  with  the  precise  determination  of  three  important 
paramagnetic  resonance  parameters  of  the  ruby  crystal  at  room  tempera- 
ture.  These  are  the  two  spectroscopic  splitting  factors  (g(.  and  gj_  ) 
and  the  crystal  field  splitting  (D)  of  the  chromium  ion  energy  levels. 

The  values  obtained  ape:   g..  -  1.98145  +  .00020,  g±=   1.98137 
+  .00067,  D  =  5.74674  +  .00080  Kmc/sec. 
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1.   Introduction 

In  order  to  determine  precisely  the  parameters  of  a  given  physical 
system,  the  experimentalist  normally  requires  a  precise  theory  describ- 
ing the  phenomena  being  observed.   For  paramagnetic  resonance  in  crys- 
tals, the  theory  which  best  predicts  the  behavior  of  the  spin  energy 
levels  of  the  paramagnetic  ions  in  an  externally  applied  magnetic  field 
is  the  quantum  mechanical  spin  Hamiltonian  developed  in  the  references 
1,  2,  5  and  7. 

The  specific  form  of  the  spin  Hamiltonian  used  here  applies  to  an 
'  ionic  crystal  field  of  medium  strength.   Medium  strength  means  that  the 
crystal  field  interaction  with  the  ion  spin  is  stronger  than  its  spin- 
orbit  coupling.   The  Hamiltonian,  used  is  also  abbreviated.   The  terms 
for  the  spin-orbit  interaction,  nuclear  spin-electron  spin  interaction, 
and  other  much  smaller  terms  have  been  omitted.   This  simplified 
Hamiltonian  is: 


^=PoH'g-S+S-D-S 


The  first  term  is  the  external  magnetic  field  interaction  with  the 
magnetic  moment  of  the  ion;  where  H  is  the  magnetic  field  strength 
vector,   g   is  the  spectroscopic  splitting  factor  tensor,  and  s   is  the 
"effective-spin"  vector  operator  which  does  not  necessarily  have  the 
same  total  eigenvalue  as  the  free  ion.   The  second  term  is  the  inter- 
action energy  of  the  crystal  field  with  valence  electrons  of  the  ion; 
where  D   is  the  electrostatic  field  strength  tensor.   The  constant   3 
is  the  Bohr  magneton. 

The  procedure  followed  here  to  find  the  energy  eigenvalues  of  this 
Hamiltonian  is  to  diagonalize  its  matrix  representation.   Therefore,  the 


development  of  the  matrix  Hamiltonian  will  now  be  roughly  outlined. 

The  quantities  in  each  term  must  first  be  referred  to  the  same  co- 
ordinate system.   The  best  choice  is  the  principal  axis  system  of  the 

A 

crystal  where  D   is  diagonal.   In  addition,  if  the  crystal  has  axial 
symmetry,  as  in  the  case  of  ruby,  the   z   axis  can  be  chosen  to  be  the 
axis  of  symmetry  and  the   x  and  y  axes  chosen  arbitrarily  in  the  per 
pendicular  plane.   In  order  to 
simplify  the  Hamiltonian  the  ex- 
ternal field  will  be  used  to 
orient  the  axes  such  that  H   lies 
in  the   x-z   plane  and  makes  an 
angle   9  with  the   z   axis.   See 
Figure  1.   The  Hamiltonian  now  be- 


Figure  1 


comes 
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The  assumption  that   g   is  also  now  diagonal  is  valid  from  group 

theoretic  arguments. 

The  convention  now  is  to  add  a  constant  energy  term,  -1/»(D   +  D 

'  3   xx     yy 

2 
+  D   )  S  (S  +  1) ,  and  regroup  the   S    terms.   This  additive  constant  has 

the  effect  of  placing  the  zero  energy  reference  level  midway  between  the 
crystal  field  splitting.   Also  since  we  have  axial  symmetry  the  follow- 
ing quantities  are  defined:   g,   =  g   =  g    ,  g„   =  g    .   After  re- 

&J.  &xx        &yy        °\\  &zz 

grouping,    the  Hamiltonian    is: 


y    =   30  H   cosB  g|1     Sz  +  PQ  H   sin6  gj_  Sx 
+D   Is]    -    1/3S(S  +   l)\     +  E(S*    -    Sy) 

where:      D   =  D        -    1/0(D        +  D      )     ,      E    ■   1/0(D        -   D      ) 
zz  /2N    xx  yy  '2V   xx  yy' 


The   constant     E      is   zero   in  axial   symmetry   since     D        ■  D 

J  J  xx    yy 

It  is  found  from  resonance  data  on  ruby  that  an  external  magnetic 
field  splits  the  paramagnetic  ion  energy  levels  into  four  levels.   There- 
fore, since  the  number  of  nondegenerate  spin  levels  must  equal  2S  +  1  , 
then  the  effective  spin  equals  3/2  •   The   S  matrices  are,  thus,  the 
4x4  spin  3/2  matrices.   That  is: 
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Closed  solutions  for  the  eigenvalues  of  the  matrix  cannot  be  obtain- 
ed except  for  the  two  cases   9=0   and   9  =  90   .   Therefore,  these 

two  sets  of  solutions  were  chosen  as  the  theoretical  equations  for  the 

* 
experimental  data. 

•    The   8=0   case  gives: 


ty(0°)= 


3/2eoHS..  +D         0  0  0 
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Which  is  already  diagonal;  therefore  the  energy  eigenvalues  are: 

El  =  3/2  Po  H  8„  +  D 
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The   9  =  90   case  gives 
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When  this  matrix  is  diagonalized  the  following  energy  eigenvalues 
are  obtained: 
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The  subscript  notation  has  been  chosen  so  that  the  0   and  90 
energies  agree  for  a  free  ion  having  D  =  0  and  g  ,  =  g..  . 

As  stated  before,  the  Hamiltonian  used  here  neglects  some  terms  and 
assumes  an  electrostatic  crystal  field.   Therefore,  the  proof  of  the  ac- 
curacy of  this  approximate  Hamiltonian  in  predicting  the  correct  energy 
levels  lies  in  the  accuracy  with  which  the  experimental  data  fits  the 
formulae.   If  the  deviations  of  the  data  from  the  equations  are  less 
than  the  experimental  errors  then  the  theory  is  experimentally  accept- 
able and  the  parameters  obtained  are  accurate  to  within  experimental 
error. 


2.   Description  of  Apparatus 

The  equipment  employed  to  make  the  paramagnetic  resonance  measure- 
ments is,  essentially,  a  research  model  of  a  conventional  paramagnetic 
resonance  spectrometer.   It  utilizes  the  microwave  balanced  bridge  tech- 
nique with  crystal  detection  and  heterodyne  amplification  of  the  reso- 
nance signal.   Figure  2  shows  a  schematic  representation  of  the  system. 

The  microwave  bridge  is  comprised  of  a  "hybrid-tee"  junction  having 
four  microwave  arms.   This  junction  divides  the  power  as  shown  by  the 
the  arrows  in  figure  2.   Two  of  the  arms  act  as  opposing  loads  and  the 
remaining  two  are  a  power  input  arm  and  a  signal  output  arm.   The  junc- 
tion divides  the  input  power  equally  between  the  loads  and  recombines 
both  reflected  signals  to  give  a  resultant  output.   One  of  the  load  arms 
is  terminated  by  an  attenuator  and  phase  adjustor  while  the  other  is 
terminated  by  a  cavity  containing  the  ruby  in  a  magnetic  field.   Conse- 
quently, if  the  phase  and  attenuation  in  the  dummy  load  arm  are  adjusted 
for  minimum  output  then  any  absorbtion  of  power  in  the  cavity  will  show 
up  as  an  unbalance  signal.   This  unbalance  is  thus  a  very  sensitive  in- 
diction  of  absorption  by  the  ruby.   Therefore,  if  the  input  frequency 
is  held  constant  and  the  magnetic  field  swept  through  resonance  condi- 
tions, an  output  resonance  signal  will  be  obtained.   This  signal  is  then 
detected  by  a  crystal,  amplified  and,  after  rectification,  displayed  on 
an  oscilloscope  for  measurement  adjustments.   Due  to  the  method  of  bridge 
balancing,  a  resonance  signal  could  be  obtained  which  had  any  character- 
istic between  pure  absorptive  or  pure  dispensive:   therefore,  a  certain 
amount  of  skill  was  required  to  obtain  the  desirable  absorption  signal. 
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Figure   2,  r 


The  microwave  power  is  supplied  to  the  bridge  by  a  reflex  klystron 
operating  in  the  frequency  range  of  12.4  to  18  kilomegacycles  per  second 
and  a  power  output  between  50  to  100  milliwatts.   The  limitation  on  the 
frequency  range  imposed  a  restriction  on  the  spread  between  data  points. 
However,  even  over  this  frequency  range,  the  different  transitions  gave 
a  spread  in  fields  strengths  which  nearly  covered  the  measurable  field 
range . 

The  heterodyne  beat  frequency,  which  is  required  for  amplification 
of  the  resonance  signal,  is  generated  by  coupling  into  the  output  arm 
of  the  bridge  a  microwave  signal  from  a  second  klystron.   The  frequency 
of  this  klystron  is  locked  onto  the  main  signal  so  that  there  is  a  dif- 
ference of  30  megacycles  per  second,  which  is  the  intermediate  frequency 
of  the  amplifiers.   This  locking  is  a  accomplished  by  using  a  30  mc  des- 
criminator  and  feeding  its  error  voltage  back  to  the  reflector  voltage 
of  the  second  klystron.   If  the  frequency  of  the  second  klystron  is 
properly  located  with  respect  to  the  main  klystron  frequency,  then  this 
feedback  will  stabilize  the  difference  frequency. 

Frequency  stability  of  the  main  klystron  is  thus  important,  since 
the  stability  of  the  entire  system  is  dependent  on  it.   It  was  found 
that  a  water  cooled  jacket  mounted  on  the  klystron  gave  very  good  stabil- 
ity (about  1  mc  drift  per  hour).   This  stability  was  good  enough  to  per- 
mit the  assumption  of  one  frequency  value  for  a  whole  series  of  field 
measurements  on  a  number  of  different  resonances. 

In  order  to  measure  the  frequency  of  the  main  signal,  a  cavity  type 
frequency  meter  is  inserted  in  the  output  arm  of  the  bridge.  This  meter 
has  a  guaranteed  accuracy  of  +  0.17o  over  a  wide  variation  in  ambiant 
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temperature  and  relative  humidity.   However,  after  some  preliminary  data 
was  reduced,  it  was  found  that  the  parameters  differed  from  values  given 
in  other  references;  and  the  accuracy  of  the  meter  was  questioned.   It 
was,  therefore,  decided  that  a  more  precise  frequency  measurement  was 
needed  and  the  organic  substance,  diphenyl  picryl  hydrazyl  (DPPH) ,  was 
chosen  as  a  natural  paramagnetic  resonance  reference.   DPPH  has  the  con- 
venient properties  of  an  isotropic  splitting  factor,  which  is  tabulated 
to  five  decimal  places,  and  an  effective  spin  of  1/2,  which  results  in 

only  one  possible  energy  transition.   The  DPPH  sample  was  about  one  bun- 
's 

dreth  of  a  gram  of  small  crystals  placed  between  a  celluloid  tape  sand- 
wich.  This  sample  was  placed  in  the  cavity  in  close  proximity  to  the 
ruby.   Use  of  this  natural  paramagnetic  resonance  reference  makes  it 
possible  for  future  investigators  to  check  or  correct  the  parameters  ob- 
tained here. 

The  magnetic  field  is  supplied  by  an  electromagnet  which  has  twelve 
inch  diameter  pole  faces  and  a  maximum  field  strength  of  about  nine  kil- 
ogauss.   The  magnet  is  mounted  on  a  rotatable  base  with  angular  gradu- 
ations, making  different  field  orientations  possible.   The  accompanying 
magnet  power  supply  has  two  stages  of  current  regulation  and,  therefore, 
has  very  good  field  stability.   The  pole  caps  of  the  magnet  have  been 
carefully  shimmed  to  give  a  high  field  homogeniety;  for  example,  at  2.5 
inches  of  radial  displacement  from  the  center,  a  3.5  kilogauss  field 
drops  0.4  gauss.   This  high  homogeniety  insured  a  negligible  broadening 
of  the  resonance  due  to  inhomogeniety  and  allowed  field  measurements  ex- 
ternal to  the  cavity  without  correction  for  the  point  of  measurement. 


The  field  strength  measurements  were  made  using  a  nuclear  spin 
resonance  gaussmeter  with  a  range  of  0.5  to  9.0  kilogauss  when  a  water 
sample  is  used.   If  an  oscillator  and  frequency  counter  are  used  in  con- 
junction with  this  gaussmeter,  an  accuracy  of  1  part  in   10   is  possible, 

As  a  final  statement  concerning  the  apparatus,  it  should  be  mention- 
ed that  the  largest  contributors  to  measurement  uncertainty  were  signal 
noise  and  bridge  balance  drift.   The  signal  noise  was  due  to  60  cycle 
pickup  and  vibrations.   Both  of  these  were  greatly  reduced  by  using  a 
direct  current  filament  supply  for  the  amplifiers  and  by  mounting  the 

-a 

entire  bridge  system  on  a  vibration  free  platform.  The  bridge  balance 
drift  was  minimized  by  constructing  the  dummy  load  arm  to  be  as  nearly 
alike,  physically,  to  the  cavity  arm  as  possible. 
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3.   Experimental  Procedure 

Before  any  measurements  could  be  made,  the  ruby  sample  had  to  be 
correctly  oriented  with  respect  to  the  magnetic  field.   The  ruby  sample 
was  in  the  form  of  1/4  inch  diameter,  1-1/8  inch  long  rod.   The  faces 
had  been  cut  very  carefully  perpendicular  to  the  crystal  axis  of  sym- 
metry and  polished  flat.   Therefore,  as  an  initial  allignment  procedure, 
one  ruby  face  was  used  as  a  reference  plane  and  its  polished  surface  as 
an  optical  mirror.   Assuming  that  the  magnet  field  was  perpendicular  to 
the  pole  faces  and  knowing  how  much  off  of  vertical  these  were,  it  was 
possible  to  set  up  a  light  beam  system  which  was  referred  to  the  hori- 
zontal.  A  cathetometer- and  lamp  were  set  up  as  in  figure  3  and  the  ruby 
was  adjusted  until  its  face  was 


pl^wc-L+o  P°!e  -Peaces 


f—? 


I  o.  flip 
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off  of  vertical  by  the  same 
amount  as  the  magnet  pole  faces. 
The  magnet  was  then  rotated  to 
a  position  where  the  ruby  face 
and  magnet  pole  face  were  approx- 
imately parallel;  this  being 
approximately  a  zero  degree 
orientation.   A  resonance  signal 
was  now  obtained  on  the  oscillo- 
scope and  the  magnet  slowly 
rotated  back  and  forth.   The 
signal  was  seen  to  move  to  the 
left  on  the  scope,  reach  a  minimum,  and  move  back  to  the  right  as  the 
magnet  was  rotated  continuously  in  one  direction.   This  minimum  was 


line 


Figue  3.   Optical  Allignment 
Setup 
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assumed  to  occur  when  the  crystal  axis  was  correctly  oriented  parallel 
to  the  field.   The  90  degree  position  could  thus  be  obtained  by  rotating 
the  magnet  exactly  90  degrees  from  this  position.   The  possible  align- 
ment error  was  less  than  0.2  degrees  if  the  assumption  was  correct  and 
the  ruby  face  had  been  accurately  cut. 

The  resonance  measurements  were  now  begun.  Three  frequencies  were 
used,  one  at  the  high  and  one  at  the  low  limits  of  the  klystron  and  one 
about  half  way  between.  At  each  frequency  the  following  transitions 
were  measured.  The  two  allowed  zero  degree  transitions,  E~  -  E-  and 
E,  -  E2  ,  and  four  90  degree  transitions,  E2  -  E~  ,  E,  -  E2  >  E2  -  E,  , 
and  E.  -  E~  .  Also  periodic  measurements  were  made  on  the  DPPH  transi- 
tion for  frequency  determination. 

The  actual  signal,  upon  which  the  measurements  were  based,  was  seen 
as  an  oscilloscope  display  of  the  absorption  curve.   This  was  accomplish- 
ed by  taking  the  detected  output  of  the  final  amplifier  and  putting  it 
into  the  vertical  deflection  input  of  the  oscilloscope  and  simultaneously 
putting  a  phase  adjustable  60  cycle  per  second  signal  into  the  horizontal 
deflection  input.   The  main  field  was  modulated  with  a  60  cycle  per  sec- 
ond field  having  a  peak  to  peak  amplitude  of  about  50  gauss.   Conse- 
quently, if  the  phase  of  the  horizontal  input  signal  on  the  oscilloscope 
was  adjusted,  a  dual  overlapping  display  of  the  resonance  curve,  as  in 
figure  4,  could  be  obtained.   If  this  absorption  curve  was  assumed  to 
be  symmetrical,  then  the  main  field  could  be  adjusted  until  the  zero  of 
the  modulating  field,  or  horizontal  oscilloscope  input,  was  at  the  cross- 
ing of  the  two  curves  and  the  main  field  would  then  be  at  the  center  of 
the  resonance.   The  symmetry  of  the  absorption  curve  was  checked  by  use 
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of  a  phase  locked  amplifier  and  a 
37  cycle  per  second  modulation 
field.   The  derivative  curve, 
thus  obtained,  was  found  to  be 
very  close  to  symmetrical  and 
therefore,  the  integral,  or  the 
absorption  curve  itself,  was 
found  to  introduce  only  about 

+  .017o  error  due  to  asymmetry.        Figure  4.   Oscilloscope  Display 
Theoretically  this  absorbtion  curve  at  room  temperature  should  be  very 
symmetrical,  but  due  to  the  method  of  bridge  balancing,  the  author  was 
not  certain  of  his  accuracy  in  obtaining  a  pure  absorbtion  signal  and 
the  derivative  plot  was  needed  for  assurance  that  only  a  very  small  dis- 
persive component  was  present. 

When  the  field  had  been  accurately  set  at  the  center  of  the  reso- 
nance, then  its  strength  was  measured  by  use  of  the  nuclear  spin  reso- 
nance gaussmeter.   This  involved  accurately  centering  the  nuclear 
resonance  signal,  as  displayed  on  an  oscilloscope,  at  the  zero  of  a  60 
cycle  per  second  modulation  field  and  measuring  the  frequency  of  the 
gaussmeter  oscillator.   This  frequency  was  determined  by  beating  a  signal 
from  a  generator  against  that  of  the  gaussmeter  and  counting  the  gener- 
ator  frequency  when  zero  beat  was  obtained.   Also  at  zero  beat  the  nuclear 
resonance  signal  was  again  examined  to  guarantee  that  there  had  been  no 
frequency  shift  caused  by  the  generator  signal. 

This  whole  measurement  procedure  was  repeated  a  number  of  times  for 
each  transition  in  order  to  obtain  an  estimate  of  the  experimental  un- 
certainty. 
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4.   Data 

In  order  to  correlate  the  resonances  measured  to  particular  energy 
level  transitions,  it  was  necessary  to  plot  the  energy  levels  in  dimen- 
sionless  coordinates.   Using  dividers  to  determine  the  positions  of  each 
transitions  for  a  constant  frequency  in  a  decending  order  of  field 
strengths  and,  knowing  approximately  what  the  absolute  scales  were,  it 
was  possible  to  identify  the  resonances.   The  energy  level  plots  are 
shown  in  figure  5. 

The  following  table  gives  the  data  after  averaging  the  repeated 
measurements  and  finding  standard  deviations.   The  values  are  the  fre- 
quencies of  the  nuclear  resonance  gaussmeter  for  the  magnetic  field  at 
the  indicated  paramagnetic  resonances.   The  ruby  used  had  a  chromium 
molar  percentage  of  0.17o.   The  ambient  temperature  was  22  ±  2   C. 


TABLE  I 
9  =  0° 


e  »  90 


Ruby  resonance 


DPPH  resonance 


E2  -E3 

El  "E2 

1/2-  -1/2 

26.9164  *  .0010  mc 
23.3063  f  .0010 
20.6576  +  .0010 

9.2632  +  .0005  mc 
5.6444  ±   .0005 
2.9840  +  -0005 

26.6087  +  .0020  mc 
23.0310  +  .0020 
20.4006  -t   .0020 

El  "E3 

E2  "E4 

1/2-  -1/2 

(not  obtained) 
4.3781  ±   .0010  mc 
2.5264  +  .0012 

15.4879  +  .0008  mc 
13.1963  +  .0005 
11.3377  t  -0015 

26.6087  +  .0020  mc 
23.0310  +  .0020 
20.4048  +  .0015 

E2  "  E3 

El  -E2 

1/2-  -1/2 

24.5077  ±  .0010  mc 
20.4065  +  .0015 
17.1962  +  .0015 

18.7616  +  .0003  mc 
15.4376  +  .0015 
13.0880  +  .0010 

26.6087  ±  .0020  mc 
23.0310  +  .0020 
20.4006  +  .0020 
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Figure   5:      Energy  Levels    in  Dirnensionless   Coordinates 
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5.   Parameter  Determination 

Taking  the  differences  of  the  energies,  as  given  in  the  introduc- 
tion, and  dividing  by  planck's  constant  (h) ,  the  desired  frequencies  of 
the  transitions,  which  were  measured,  are  obtained. 

For  zero  degrees: 

B  B 

V-3  '  "h  H  g"   ;   "l-2  "  HI  H  «M  +  2D 

Where  D   is  now  in  units  of  frequency 

For  ninety  degrees: 
8 


For  DPPH:  « 


»-  -THg 


The  constants  used  are: 

3 


1  moa  _i_   nnnn/   mc/sec  (American  Institude  of 

—r     =      1.39966  +   .00004  -  „.    „   ,,   ,  N 

h  -  gauss   Physics  Handbook) 
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g  =  2.00354  +   .00003  (Reference  6) 

Also  since  the  field  was  measured  in  terms  of  the  proton  resonance  fre- 
quency in  water,  the  following  was  used: 


H(gauss)  -  (234.861  +  .010)  V      (mc/sec) 

Inserting  the  constants  and  using  only  proton  frequencies  for  DPPH 

fields  and  ruby  fields,  for  the  same  microwave  frequency,  the  zero  degree 

equations  are  solved  for  g..  and  D  . 

V   (DPPH) 
for   (E2  -  E3):   g|(   =  (2! 00354  +  .00003)   /  (fuby) 

for   (Ex  -  E2):   D  ■  (329.308  +  .020)   (u  (DPPH) —  up  (ruby)] 

For  the  90   frequency  equations  the  following  substitutions  are  made: 

3  H 
A  ■  -£ —  =  (328.726  +  .020)  up  (ruby) 

B  =     V  =   (658.615  +  .040)  u  (DPPH) 

and   the   equations   were   solved   for     D     as   a   function  of     g.    ,   A  and  B. 


(B  +  Ag,  )  /  B(2Ag_L  +  B)    -    3(Ag,)2 

for        (E,    -  EQ):      D   = 


2        "3'*  2  W  B(2Agx  +  B) 


for        (Ex   -  E3):      D   =     \  ./  B~     ■   e(Agj_  )"         ■  Agx 


for        (E2    -  E4) 


:      D   =     |  ^B2    -    3(Agx  )2        +  Agx  j 


(B        AgjL  )           /B(2Ag±  -    B)    +   3(Agj_  )2 
d1       E2):      D<      — \     B(2Ag      -   B) 
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a)   Zero  Degree  Solutions 

Inserting  the  data  for  the  E„  -  E„   transitions  into  the  g.. 
equation,  the  following  values  were  obtained. 

g(t  =  1.9806  +  .0002   for   v         =  17.525  Kmc/sec 

g„  -  1.9799  +  .0002   for   i?2_3  =  15.169  Kmc/sec 

gu  -  1.9786  +.0002   for   u    =  13.436  Kmc/sec 

These  values  do  not  agree  within  experimental  error  and  there  is  an 
evident  frequency  correlation.   This  frequency  correlation  suggests  a 
misallignment  of  the  ruby  axis  and  the  magnetic  field.   In  order  to  esti- 
mate the  amount  of  misallignment,  data  from  a  previous  run,  in  which  a 
frequency  meter  was  used,  is  plotted.  (See  figure  6)   This  data  is  from 
measurements  of  an  E„  -  E„   transition,  at  13.597  Kmc/sec  as  a  function 
of  angle  about  zero  degrees.   The  crystal  was  a  different  one  from  that 
used  here  and  had  a  chromium  percentage  of  .  057o,  therefore,  the  misallign- 
ment encountered  here  is  a  cutting  error  of  the  0.17»  Cr  crystal  face. 
The  data  has  been  reduced  in  the  form  of  g  values,  obtained  using  the 
zero  degree  equation,  versus  angle  from  zero. 

0°     g     1.9811  +  .0010 

+   1°  gx0   ■   1.9790  +   .0010 

+  2°  g2°   =   1.9727  +   .0010 

+  3°  *u°   =   1.9643  +   .0010 

o  ±|w         ±Z"        ±3"        xi" 

+  4°  g,°  =   1.9517  +   .0010 


Figure  6.   g  versus  angle 
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The   accuracy  of   this    data    is   not   as   good  as    the   present   data   but    is 
still   useful.      The   curve   closely   fits   a  parabola  having   the    formula: 

g0  *  g,|      "    (-0020)    92   =   1.9811   -    ,.0020)    92 

Where   9   is  in  degrees . 

Comparing  the  value  obtained  at  13.436  Kmc/sec  to  this  gives  an 
angle  of:   9  =-  1.25°  . 

Rather  than  attempt  to  allign  the  ruby  more  accurately,  another 
approach  was  used.   Since  the  data  obtained  is  accurate  for  some  small 
angle  approximately  equal  to  one  degree,  the  energies  of  the  zero  degree 
Hamiltonian  are  modified  to  take  into  account  a  small  angle  from  zero. 

The  approach  is  to  take  the  general  Hamiltonian,  substitute   9   for 

2 
sin9  and  (1-9  /2)  for  cos9  and  use  perturbation  theory  on  the  resultant 

Hamiltonian.   See  Appendix  I  for  the  derivation  of  these  small  angle 

energies.   The  equations  for   g   and  D  are  put  in  an  approximate 

form  in  order  to  use  iterative  techniques.   It  is  assumed  that 

g,   =  g,j  -  gn   in  the  correction  terms. 


hv> 


2-3 


+  l  So 


2  D  h 

PoSoH 


D  =  x 


1-2 


"  Si 


H  +  2  go  T  H 


1  - 


1  + 


2  D  10 

eos0Hy 


From  the  previous  data  a   g   =  1.9813  was  obtained;  therefore,  this  is 
used  for  the   g   value.   Also  if  the  present  data  is  used  in  the  zero 
degree  equations  for  D  ,  the  values  obtained  do  not  show  any  misallign- 
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merit  effect  and  an  average  value  D  =  5.7462  Kmc/sec   is  obtained.   This 
is  used  in  the  small  angle  equations  in  the  correction  term. 
Using  the  E~  -  E„   data,  the   g  .   equations  are: 

(1)  g)(  =  (1.98064  +  .00020)  +  (2.2397)  92  for  V         =  17.525  Kmc/sec 

(2)  gj(  =  (1.97988  +  .00020)  +  (3.9914)  92  for  1>2   =  15.169  Kmc/sec 

(3)  g(l  =  (1.97861  +  .00020)  +  (8.0195)  92  for  u£_3  =  13.436  Kmc/sec 

2 
At  this  point  it  is  interesting  to  note  that  the  constant  8.0195/rad 

2 
in  equation  (3)  is  . 00243/d§g  and  compares  nicely  with  the  experimen- 

2 
tally  obtained  value  of  .0020/deg  ,  which  was  for   uo  o  =  13.597  Kmc/sec. 

Now,  there  are  three  equations  and  two  unknowns   g.,  and   9   ; 

therefore,  three  values  of  g..  can  be  obtained  for  an  averaged  value 

of   6  .   The  average  value  of   9   is: 

92  =  .0003668  +  .0000500  radians2 

or        9  =  1.097  +  .080  degrees 

from   (1)  g()  =  1.98146  +  .00020 

(2)  g|)  =  1.98134  +  .00020 

(3)  g(|  =  1.98155  +  .00020 

The  average  of  these  three  is: 


g,   =  1.98145  +  .00020 


Using  this  value  of  g    in  the  D  equations,  the  same  process 
can  be  usede 

(4)  D  =  (5.74565  +  .00080)  +  (2  9673)  92 

(5)  D  =  (5.74605  +  .00080)  +  (2.1677)  92 
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(6)   D  =  (5.74625  +  .00080)  +  (1.2469)  92 

The  average  angle  is: 

92  =  .0003554  +  .0001300  rad2 

6  =1.080+0.20  degrees 

from   (4)   D  =  5.74670  +  .00080  Kmc/sec 

(5)  D  =  5,74682  +  ,00080  Kmc/sec 

(6)  D  =  5.74669  +  .00080  Kmc/sec 

The  average  of  these  three  is: 


D  =  5.74674  +  .00080  Kmc/sec 


b)   Ninety  Degree  Solutions 

The  determination  of  g  ,  and  D   from  the  ninety  degree  equations 
is  more  involved  than  the  simple  solution  of  a  set  of  linear  algebraic 
equations  like  the  zero  degree  equations.   The  ninety  degree  equations 
are  a  set  of  quadratic  algebraic  equations.   Using  the  available  data 
there  are  eleven  equations  and  the  two  unknown  parameters.   One  could 
assume  the  same  D  value  as  already  obtained  from  the  zero  degree  data 
but  this  would  not  allow  an  independent  check  of  the  accuracy  of  the 
ninety  degree  equations.   Therefore,  the  equations  were  solved  simulta- 
neously using  the  following  method.   If  eleven  simultaneous  equations 
in  two  unknowns  are  solved,  fifty  five  intersection  values  are  obtained. 
Since  this  is  a  lengthy  procedure  the  D  equations  were  programmed  in 
an  IBM  1604  computer  and  the  values  of  D   for  each  of  the  eleven  equa- 
tions at  values  of  g    ,  running  from  1.970  to  2.000  in  steps  of  .001, 
were  determined.   An  index  i  was  used  to  label  the  data  sets  and   j   the 
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values  of  g.   .   The  values   i  =  1,2,3  are  for  E..  -  E?   data,  i  =  4,5,6 

for  E2  -  E„   data,  i  =  7,8,9   for  E2  -  E,   data,  and  i  =  10,11,12  for 

E,  -  E„  data.   Where  each  triplet  is  started  with  the  low  value  of  B. 
13  i 

and  increased  with   i  .   The  computer  was  instructed  to  label  each  column 

of  D.  .   values  as   D(i)  .   Table  II  gives  the  program  print  out.   There 

was  no  data  for   i  =  10  and  the  zeros  were  printed  if  the  value  of  D.  . 

ij 

was  imaginary.   If  this  data  is  plotted  as  in  figure  7,  the  scattering 
of  the  intersection  points  can  be  seen.   The  results  are  better  than 
this  plot  suggests  since  each  intersection  does  not  carry  the  same  weight 
in  determination  of  g,   and  D  average   This  is  because  data  error 
gives  the  lines  a  certain  width  and  thus  the  intersection  coordinates 
are  affected  not  only  by  this  width  but  by  the  angle  between  the  lines 
at  intersection  and  their  orientation  with  respect  to  the   g  .  and  D 
coordinate  axes.   Appendix  II  gives  the  derivation  of  the  weighting  factor 
formulae  which  were  used  in  the  averaging.   The  values  of  D  and  g  , 
at  each  intersection  were  determined  by  locating  the  interval  between 
which  the  two  curves  intersect,  in  Table  II,  and  using  straight  line  pro- 
portionality between  the  four  D  points. 

The  final  results  are  given  in  the  next  section. 
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G  01  D2  C3  CU  D5         06 

1.970000  5.8U5222  5.837711  5.8495C4  5.898652  5.847036   5,822137 

1.971000  5.856665  5.829438  5.8U10C4  5.885*18  5.838501   5.816448 

1.972000  5.828119  5.821181  5.832527  5.872942  5.8299U3   5.810746 

1.973000  5.819585  5.812939  5.824071  5.860023  5.821 363   5.805C32 

1.974000  5.811062  5.8C4713  5.815637  5.847060  5.8127c0   5.7993C5 

1.975000  5.802550  5.796502  5.807225  5.834053  5.8C4134   5.793566 

1.976000  5.79405C  5.788307  5.798fi33  5.821003  5.795485   5.787815    i 

1.977000  5.785561  5.7-80127  5.790463  5.8C790R  5.786813   5.782049 

1.978000  5.777083  5.771962  5.782114  5.794768  5.778118   5.776271 

1.979000  5.768617  5.763813  5.773786  5.781584  5.769599   5.770481 

1.98GC00  5.760161  5.755678  5.765478  5.768354  5.760657   5.764678    '; 

1.981000  5.751717  5.747559  5.757^91  5.755079  5.751892   5.758862 

1.982000  &.7U32a|  1.739454  5,74g925  5.741757  5.745103   5.7;>iCJi 

1.985000  5.734860  5.731364  5.740679  5.72839C  5.734290   5.7u?19f 

1.98U000  5.726449  5.723289  5.732453  5.714976  5.725^54   5.741337 

1.985000  5.718048  5.715228  5.7242^7  5.701515  5.716593   5.735469 

1.986000  5.709657  5.7C7182  5.7160cl  5.688007  5.7C7708   5.729588 

1.987000  5.701278  5.699150  5.707895  5.674451  5.698800   5.723695   i' : 

1.988000  5.692909  5.691133  5.699748  5.66C847  5.689866   5.717783 

1.98900C  5.684550  5.683130  5.691621  5.647196  5.680909   5.711868 

1.990000  5.676202  5.675141  5.683514  5.633495  5.671927   5.7C5934    / 

1.991C00  5.667865  5.667166  5.675426  5.619745  5.662920   5.699988 

1.992000  5.659538  5.659205  5.667357  5.6C5947  5.653889   5.694028 

1.993000  5.651221  5.651258-  5.6593C7  5.592098  5.644832   5.688C55 

1.994000  5.642915  5.643325  5.651276  5.578199  5.635751   5.662C68 

1.995000  5.634619  5.635406  5.643264  5.56425C  5.62664U   5.676063 

1.996000  5.626333  5.627501  5.635271  5.550250  5.617512   5. 670C54 

1.997000  5.618057  5.6196C9  5.627296  5.536198  5.608355   5.664027 

1.998000  5.609792  5.611731  5.61934C  5.522095  5.599172   5.657986 

1.999000  5.601536  5.603866  5.6114C2  5.50794C  5.589963   5.651932 

G  D7  C8  D9  D11  012 

1.970000  6.169410  5.931U98  5.84385C  5.758443  5.750374 

1.971000  6.138282  5.916813  5.836245  5.757297  5.749804 

1.972000  6.106093  5.901945  5.828593  5.756151  5.749233 

1.973000  6.072736  5.886888  5.820894  5.755004  5.748663 

1.974000  6.038090  5.871634  5.813U8  5.753857  5.748092 

1.975000  6.002004  5.856179  5.805353  5.752710  5.747521 

1.976000  5.964501  5.840515  5.7975C9  5.751562  5.746950 

1.977000  5.924763  5.824635  5.789615  5.75C415  5.746380 

1.978000  5.883115  5.3C8531  5.781671  5.749267  5.745*09 

1.979000  5.859009  5.792195  5.773676  5.748119  5.745237 

1.980000  5.791987  5.775619  5.765629  5.746970  5.744666 

1.981000  5.741429  5.758792  5.75753C  5.745822  5.744095 

1.982000  5.686^(56  5.741706  5.749577  5.744673  5.743524 

1.983000  5.625751  5.724349  5.741170  5.7u3524  5.742952 

1.984000  5.557172  5.706710  5.732908  5.7U2375  5.742381 

1.985000  5.476772  5.688778  5.724590  5.741225  5.741809 

1.986000  5.375481  5.670538  5.716216  5.74C075  5.741237 

1.98700C  5.216465  5.651978  5.707784  5.738925  5.7^0666 

1.988000  .000000  5.633082  5.699294  5.737775  5.740094 

1.989C00  .000000  5.613834  5.690744  5.736624  5.739522 

1.990000  .000000  5.594216  5.682135  5.735474  5.738950 

1.991000  .000000  5.574209  5.673464  5.734322  5.738378 

1.992000  .000000  5.553791  5.664731  5.733171  5.737806 

1.993000  .000000  5.532941  5.655935  5.732020  5.737233 

1.994000  .000000  5.511632  5.6U7C74  5.730868  5.736661 

1.995000  .000000  5.489837  5.638149  5.729716  5.736089 

1.996000  .000000  5.467525  5.629157  5.728564  5.735516 

1.997000  .000000  5.444662  5.620097  5.727411  5.734944                                   M 

1.998000  .000000  5.421210  5.61C969  5.726258  5.734371                                   ft 

1.999000  .OOOOOC  5.397127  5.601770  5.725105  5.733798 

TIME,       0    NINUTES    AND  14    SECONCS 

Table  II:   Data  Print  out  for  Ninety  Degree  Solutions 
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^Figure    7:      Plot  of  Ninety  Degree  Computer  Data 


24 


■         :5:79.,Vi)V;5,,80, 

i     ',"■.•..■  Wi    ■ 
1    ;        ■'       '    ■  ■   i III'  , 

/  -  •■■■■■  .:a -.',—/  '     ■ 


6.   Results 

From  the  ninety  degree  data: 


gx  =  1.98137  +  .00067 

D   =  5.7489  +  .0061  Kmc/sec 


From  the  zero  degree  data; 


g   =  1.98145  +  .00020 

D   =  5. 74674  +  .00080  Kmc/sec 


The  D  values  agree  and  the  g  values  are  equal  within  experimen- 
tal error. 

These  g  values  differ  quite  significantly  from  those  obtained  by 
A.  A.  Manenkov  and  A.M.  Prokhorov  (see  reference  4)  which  were 

g(|  =  1.9840  +  .0006,  gx  =  1.9867  +  .0006  and  D  =  5.732  Kmc/sec. 

Their  report  did  not  explain  in  detail  how  these  were  obtained  except 
that  they  used  the  same  Hamiltonian  as  used  here.   Also  another  author, 
J.  E.  Geusic  (reference  3),  measured  these  parameters  and  gives: 

S ||  =  2.003  +  .006,   gj_  -  2.00  +  0.02   and  D  =  5.78  +  .03  Kmc/sec 

which  disagrees  with  the  g..  and  D  of  Manenkov  and  Prokhorov  and  also 
the  g   obtained  here.   Geusic  used  essentially  the  same  Hamiltonian 
and  similarly  measured  his  data  at  zero  and  ninety  degrees.   However  he 
stated  that  his  agreement  with  calculated  values  at  90  was  only  good 
to  1%  . 

In  order  to  check  the  accuracy  with  which  the  calculated  values  of 
the  resonance  frequencies  fits  the  experimental  values,  the  following 
tables  were  compiled  using  the   g  values  obtained  and  the  best  D  value, 
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(All  frequencies  are  in  Kilomegacycles  per  second) 
8  =  1.10  +  .08° 


transition 

(calculated) 

(measured) 

2-3 

17.525  +  .002 

17.525  +  .002 

15.169  +  .002 

15.169  +  .002 

13.435  +  .003 

13.436  +  .002 

1=2 

17.525  +  .002 

17.525  +  .002 

15.168  +  .002 

15.169  +  .002 

13.436  +  .002 

13.436  +  .002 

e  =  90 


0 


transition 

(calculated) 

(measured) 

4-3 

17.525  +  .001 

17.525  +  .002 

15.171  +  .001 

15.169  +  .002 

13.432  +  .001 

13.436  +  .002 

3=2 

17.523  +  .001 

17.525  +  .002 

15.168  +  .001 

15.169  +  .002 

13.428  +  .002 

13.436  +  .002 

3-1 

17.529  +  .002 

17.525  +  .002 

15.167  +  .001 

15.169  +  .002 

13.434  +  .002 

13.439  +  .002 

4-2 

15. ,172  +  .002 

15.169  +  .002 

13.445  +  .002 

13.439  +  .002 

The  agreement  is  very  good  for  the  small  angle  values,  but  the  ninety 
degree  values  show  a  maximum  disagreement  of  067„  which  is  slightly  out- 
side of  the  estimated  experimental  error.   This  disagreement  could 
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possibly  be  caused  by  a  misallignment  at  ninety  degrees  or  an  actual 
theoretical  disagreement.   An  approximate  perturbation  Hamiltonian  can 
be  obtained  for  a  small  angle  from  90  but  the  calculations  are  much 
more  tedious  than  for  the  zero  degree  case.   Regardless  of  the  source  of 
the  disagreement,  the  g,  value  obtained  by  a  more  accurate  method  would 
probably  not  vary  from  the  value  obtained  here  by  more  than  approximately 
0.1%  . 
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APPENDIX   I 
Small  Angle  Hamiltonian  and  Energies 
2 


%    - 
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•|PoS„H(l-e/2)+D 


Where  the  following  substitutions  were  made: 


cos 


e  =-i  -  ez/, 


in  9  -   9 


sin 


Splitting  into  an  unperturbed  and  a  perturbed  Hamiltonian 


tf-Jf„  +#' 


W- 


fe0g„H+D 


I      o 


2PoS„H-D 


lPoS„H92 

^ogj.ne 

i»oS.i«e 

-  i0og„He2 

0 

30gjH6 

0 

0 

0 

0 

o    H- 

0 

0 

Pc 

gj.H9 

X 

g„He2 

2^o 

e  H9 

0 

0 

0 

■  h 

3g..H+I) 

0 

0 

r¥ 

d8j.hG 

¥ 

^«HQ2i 
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The  zero  angle  unperturbed  energies  are: 
Eol   =   3/2  Po  SllH  +  D 
Eo2   =  Y'l   3o  8„H  -  D 
Eo3  -  -1/2  ^o  8„  H  -  D 
Eo4  =  "3/2  Po  S„  H  +  D 

The  perturbed  energies  are  obtained  using: 

4  /  / 

TTlY.fl  k£ 


K     -7fvw  + 


k         '  '  kk        /  E  ,      -  E .  » 

j  ok  o* 

4=1 

3   a  .2  3/4(Po  8j.H   Q)2 


Ex      =   -   4  3o  8„H   6     +       3q  g     h  +  2D 


"  4  Mo  H  )  [  —   -    eog(iH  +  2D 


,  V.  I|02    ,    3/4^°  »1H  fl)2      ,      (eo^H9)2 

E2     --4Po8)(H9     +  _$Q         R  _   2D         -f         po  -  H 


jtjl  3^o  §x  H  4§.L 

•   4    (gl  Po  H)      [    gj_     +     30  g     H  +  2D      '        g„ 
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e,.     -   x  3n  g„H  ez   + 


'4      4  Ko  &  ||  -PQ  g|(H  +  2D 

"  4  <3o  §1  H)   I  —  "   3og||H-  2D  '   9 
The  energies  are,  thus,  to  second  order  approximation 
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APPENDIX   II 
Weighting  Factors  for  Ninety  Degree  Data 

The  following  weighting  factor  equations  are  for  an  intersection  of 
two  straight  lines  which  have  the  same  thickness  representing  experimen- 
tal uncertainty.   Therefore,  the  two  assumptions  which  are  made  concern- 
ing the  actual  ninety  degree,   g   versus  D  ,  curves  are  that  over  the 
interval  between  known  points  the  curves  are  straight  lines  and  that  all 
the  curves  have  the  same  experimental  error  thickness.   The  first  assump- 
tion is  very  good  since  the* curves  are  nearly  straight  over  30  times  the 
smallest  interval,  and  the  second  assumption  was  verified  to  be  satis- 
factory by  computing  the  actual  uncertainties  in  g,   and  D   for  all  the 
curves.   One  factor  remains  arbitrary,  which  is  the  relative  sizes  of 
the  D  and  g,   scales  to  each  other.   In  order  to  get  the  most  variation 
among  the  weighting  factors  and  to  use  scales  which  make  the  error  bars 
for  g  ,  and  D  about  the  same  length,  the  D   scale  in  Kmc/sec  has  to 
be  l/10th  the  size  of  the  gj_  scale. 

Based  on  the  foregoing  statements  the  following  is  developed.   Fig- 
ure 8  shows  the  geometrical  relationships. 
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D 


Figure  8.   Intersection  of  Two  data  lines  having  unit  width 


From  figure  8: 


Weighting  factor  for  intersection  value  of  g  .   equals 


kgJ/     cos (3  -  oy2) 


And   for   intersection  value   of     Dt 


W(D) 


smcc 


in(3  +  a>/2) 


In  terms  of  the  slopes  of  the  two  lines 


B   aj2      =   tan' 


1  /AD 


^g 


1/  a 
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•♦•/.-  -Mil) 


lL/b 


a     =      tan 


1  /AD   \  .      -1     /AD    \ 

ta)„  -tan  tal 


The  Ag,  and  AD  values  are  obtained  from  the  data  in  table  II 
(Agi  =  .001) 
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